The Single-Electron-Box and the Helicity Modulus of an inverse square XY-Model 
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We calculate the average number ol electrons on a metal- 
lic single-electron-box as a lunction ol the gate voltage lor 
arbitrary values ol the tunneling conductance. In the vicin- 
ity ol the plateaus the problem is equivalent to calculating 
the helicity modulus ol a classical inverse square XY-model 
in one dimension. By a combination ol perturbation theory, 
a two-loop renormalization group calculation and a Monte- 
Carlo simulation in the intermediate regime we provide a com- 
plete description ol the smearing ol the Coulomb staircase at 
zero temperature with increasing conductance. 
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In recent years the field of single electronics, which is 
based on the suppression of tunneling in small metallic 
islands or quantum dots by the Coulomb blockade, has 
gained enormous interest Jl|,§l. The most elementary de- 
vice where these effects can be studied is the so called 
single-electron-box, first realized by Lafarge et al. ||. It 
consists of a small metallic island which is connected to 
an outside lead by a tunnel junction and is coupled capac- 
itively to a gate voltage V. By elementary electrostatics 
the classical Coulomb energy for a given integer number 
n of additional electrons on the box is 101 



U 



(n - n x y 



c+c Cj 



(1) 



is an ef- 



up to an irrelevant constant. Here U = 
fective single electron charging energy and n x = C g V/e 
the continuous polarization charge induced by the gate. 
Obviously, on varying n x , the actual integer value of n is 
the one minimizing E n . As a result, n(n x ) is a staircase 
function with unit jumps at n x = | (mod 1). The basic 
requirements for such a device to work are twofold: First, 
it is obvious that the temperature T (we set fcs = 1) has 
to be much smaller than the relevant charging energy. 
This point is easily taken into account by considering a 
thermal distribution of the energies E n || . The thermal 
average (n) (n x ) will then approach the simple straight 
line atT> [/ (note that n is measured over a sufficiently 
long time interval, giving a continuous (n) even though 
the number of additional electrons in the box is an integer 
at any given instant of time). Secondly, however, there 
is an intrinsic broadening of the staircase even at T = 
since the tunneling probability through the junction is 
necessarily finite. As a consequence the number of elec- 
trons in the box is not strictly conserved and the variable 
n exhibits quantum fluctuations, which are neglected in a 



simple electrostatic description. A quantitative measure 
of these fluctuations is provided by the average tunneling 
probability at the Fermi energy e p which determines the 
dimensionless tunneling conductance M 



.9 = K \t\e F PboxPlead 



h 



Ae 2 R t 



(2) 



Here t is the transfer matrix element (see (||)) and p are 
the densities of states at Sf- Obviously the perfect stair- 
case at T = is only realized in the limit g — ► 0, while 
we expect that (n) (n x ) — > n x even at zero temperature 
if g 3> 1. In the following we will give a quantitative 
description of the behaviour of (n)(n x ) near the plateaus 
at n x = (mod 1) for arbitrary values of g. It is shown 
that the slope x = at = is directly related to 
the helicity modulus of a classical 1-d inverse square XY- 
model. Combining perturbation theory to order g 2 , a two 
loop RNG calculation for g ^> 1 and a MC simulation in 
the intermediate regime, we are able to fully describe the 
crossover from perfect plateaus x(<7 —> 0) — > to the com- 
plete absence of the Coulomb blockade x(s ^ 1) — * 1- 
In particular the RNG result gives a definitive solution 
to the long-standing controversy over the preexponential 
factor at strong coupling. 

We start from a tunneling Hamiltonian 0] 



U 



H = —(h - n x f + tabclcb + h.c, + H 



(3) 



describing the Coulomb energy of the box and the trans- 
fer of electrons between states b and a from the box (in- 
dex b) to the lead (index a) and vice versa. The contri- 
bution Hq is the Hamiltonian of noninteracting Fermions 
on both sides of the junction which act as reservoirs. The 
Coulomb interaction is only incorporated by the classical 
capacitive energy, which is a good description for metallic 
systems pLEa. In the following we shall employ an effec- 
tive model for the thermodynamics of the box which is 
obtained by integrating out the fermionic degrees of free- 
dom. Using a second order cumulant expansion in the 
transfer term, which is appropriate in the experimentally 
relevant limit of a large number of conductance chan- 
nels, the reduced partition function Z = Tr cxp (~^ g ) can 

in, z « 
be written as a path integral M 



Z(n x ) = 



d6 



V6expl-S[9] 




(4) 
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over a compact angular variable 6 conjugate to the in- 
teger n. Here a; is a dimensionless coordinate of a 1-d 
system with length L — f3U where L — > oo as the physi- 
cal temperature approaches zero. The action is given by 
(-L/2 <x< L/2) 



S[9] 



dx 
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dxdx 



„: n 2 ( e(x)-e(x') , 

(x — x 1 ) 2 



(5) 



The long range part of the interaction has already been 
written in the form appropriate for L — > oo. Introducing 
a two component unit spin S(x) = (cos(0(:r)), sin(#(a:))) 
at any point of the line, the action S[6] is just the classical 
energy of an XY-model with an inverse square interac- 
tion proportional to the conductance g. The first term, 
which arises from the classical charging energy ^V 2 is 
then identical with the spin wave approximation to a 
short range interaction. Finally, the external charge n x 
acts like a purely imaginary external torque on the XY- 
model. Since J dO = 1-nm determines an integer winding 
number m G Z which is a topological invariant for each 
configuration 6(x), the external charge appears as a pure 
boundary term. Defining the free energy per length by 
f{p>x) = —lnZ(rix)/L the average number of electrons 
in the box can be expressed as 



(n) 



df{n x ) 
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The fact that n x only arises as a phase factor e 27rim ™ x 
shows that Z(n x + 1) = Z(n x ) quite generally. Thus, all 
quantities are periodic in n x with period one which allows 



to restrict the discussion to the intervall — \ 



<n x < 



In 



the following we confine ourselves to the zero temperature 
(i.e. thermodynamic) limit L — > oo and to the vicinity 
of n x = 0. Taking to be a real torque for 

the moment, a finite value of m x will induce a nonzero 
average gradient of the phase 



lim 



1 
L 



dx — ) r 

dx 
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(7) 



which is linear in m x in the limit m x — > 0. The associated 
torsional rigidity 7 is then precisely identical with the 
helicity modulus as defined by Fisher et al. S. It may 
be obtained from f(n x ) via 



1 

7 



d 2 f{n x ) 



n x =0 



(8) 



which is a measure for the sensitivity to a change in the 
boundary conditions. Using (||) the slope of the Coulomb 
staircase near n x = is related to the helicity modulus 
by 



X 



d(h) 



dn x 
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(9) 



In the trivial limiting case 3 = this describes the ex- 
pected result = 0j i- e - perfect plateaus. Indeed if 
g = 0, the helicity modulus is equal to one, being just 
the coefficient in front of the 5 term [||. In order to 

describe the behaviour at finite g we apply three different 
methods. 

Perturbation Theory: While the nonlinearity of the 
sin 2 in the action (JsJ) makes an exact evaluation of the 
path integral impossible, we may expand the long range 
contribution down to second order in g. Evaluating the 
resulting averages (exp(±i0)} with the remaining Gaus- 
sian action, the free energy can be calculated up to order 
g 2 . After a straightforward but tedious calculation we 
obtain 



x(g) = c x g + c 2 g 2 + 



(10) 



with ci = \ and C2 ~ —0.052. For the coefficient c-i we 
have evaluated a remaining definite integral numerically. 
The result ( |l7j| ) is identical with a previous calculation by 
Grabert who has used direct perturbation theory to 
fourth order in t in the original Fermionic Hamiltonian 
(Q). The agreement to order g 2 which we have verified 
to 8 digits, confirms that the inverse square XY-model 
employed here is a correct representation for the reduced 
thermodynamics of the original model. 

Renormalization Group: To obtain the exact be- 
haviour of the helicity modulus at large values of 9, 
we use the RNG. Indeed the limit g ^> 1 has previ- 
ously been treated by approximate instanton calculations 
[|7],||. However, they give different results for the pre- 
cxponential factor of the effective charging energy which 
is essentially the inverse of the helicity modulus. Here 
we will show that a definitive solution of this problem 
may be obtained by a two loop RNG, which uniquely de- 
termines both the exponent and the g-dependence of the 
prefactor of the correlation length £ in the limit g 1. 
Indeed the helicity modulus is directly proportional to 
the correlation length. To see this we define 

Pm = 2nZ(n x = 0)- 1 V9exp{-S[0]} (11) 

Jo 

as the probability for a given winding number to. It 
is then straightforward to show that the inverse helicity 
modulus 
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measures the normalized variance of the winding num- 
ber with respect to the probability distribution (pi]). At 
j > 1 a given value of (to 2 ) requires a system size which 
is larger by a factor of £(g) than that at g of order one 
where £(1) ~ 0(1). Therefore, by applying ( |l2|) we have 
7(5) ~ £(#)■ I n order to determine £(<?), we use the fact 
that d = 1 is the lower critical dimension of the inverse 
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square XY-model [||l0|. Since the kinetic energy term 
in (||) is irrelevant at g 3> 1 and T e ff — n 2 /g is the ef- 
fective temperature of our classical XY-model, we may 
perform a d = 1 + e - expansion around an ordered state 
at e > 0, which is effectively a low temperature expan- 
sion. It is convenient to generalize the XY-spin S(x) to 
a 0{n) spin S parametrized by pJ| 



8(X) = (n(x),y/l-TP{xj) 



(13) 



Here the Hi(x), i = 1, ...n— 1 are Goldstone modes whose 
expectation scales like (II 2 ) ~ T e ff. We thus expand the 
long range part of the action in powers of II. In Fourier 
space and with H as an external magnetic field which 
serves to regularize the infrared divergencies, the action 
takes the form (d = 1) 



dq- 



H 
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u q u q 
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up to terms of order II 8 . Following the method of Amit 
jllj we calculate the /3-function by field theoretic renor- 
malization using dimensional regularization. The two 
point function at zero external momentum is given by 
(e - 0) 



T (2 \0,H) 



H 



T, 



eff 
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+T eff H 1+2 * 



(n-l)(n-2) 3(n-l) : 



4e 



8f 2 



in two loop order. It can be made finite by introducing 
renormalized parameters [fLOl 



t = K £ Z- 1 T eff , h = Z^H 



and fields 



u R = z in , r ( w 2) = zr (2) 



The renormalization constant Z turns out to be 



(16) 
(17) 



Z = 1 



n-1 



-t + 



(n-1) 2 



n-1 



2c 



t 2 + 0(t 3 ) (18) 



The resulting /3-function for the renormalized tempera- 
ture is then given by 



(3(t) =et - (n - l)t 2 - (n - l)t 3 + 0(t 4 ) 



(19) 



This implies that under a reduction A — » Acxp(— I) of 
the cutoff the parameter g _1 ~ T e // at e = scales like 



dg- 1 
dl 



(20) 



By integrating this differential equation, we find that the 
associated correlation length diverges like 



> 1) =c(.g) 5 - 1 exp(2 ff ) 



(21) 



with a function c{g) which is finite as g — -> oo. The ex- 
ponential behaviour is typical for a system at its lower 
critical dimensionality, and is also obtained in an instan- 
ton approach However, the prefactor proportional 

to g^ 1 which is fixed by the coefficient of the two loop 
contribution in (pc|) and which implies that 

(22) 



X(9 > !) = 1 - cgexp(~2g) 



is quite different from the g 2 M or g 3 g] prediction of 
the instanton calculation. In fact, a very similar situa- 
tion arises in the closely related 0(n) nonlinear er-model 
in two dimensions. Due to the scale invariance of the ac- 
tion there are instantons of arbitrary size and the calcu- 
lation cannot be controlled in the thermodynamic limit. 
It is only the two loop RNG which allows to determine 
the correct prefactor of £ fl2|| , although - in contrast 
to the instanton results - it does not fix the numerical 
constant c. It should also be pointed out that the finite 
correlation length ( ^l| ) does not imply an exponential de- 
cay of (S(a;)S(0)). Indeed it can be shown |l3|] that this 
correlation function asymptotically decays like 1/x 2 for 
all < g < 7r 2 /2 and - possibly - even more slowly 
for larger values of g. For our discussion of the helicity 
modulus, however, the detailed behaviour of (S(x)S(O)) 
is irrelevant. 

Monte Carlo simulation: In order to bridge the gap 
between the perturbative result ( |l0| ) valid for small g 
and the asymptotic behaviour ( p2[ ) we have performed a 
Monte Carlo simulation of the inverse square XY model 
(0) including the short range interaction term for values 
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15| 



of g between 0.1 and 5. Following previous work 
we have sampled the winding number probabilities p m 
defined in (|ll| ) using the standard Metropolis algorithm 
with periodic boundary conditions on a discrete chain 
with up to 2000 spins. We have checked carefully that 
further increase in the system length does not change 
our results for the helicity modulus. To estimate the sta- 
tistical error we have performed about 40 runs for each 
choice of the paramters g and L. The calculations were 
done with HP 700 workstations and a CRAY T90 and 
took about one hour of CPU time per run. The numer- 
ical data are shown in Fig. [l]. Evidently the slope of 
(n) versus n x follows the perturbative result (10) closely 
up to values around g ~ 1 and finally approaches x = 1 
exponentially fast as predicted by (22). Assuming that 
the asymptotic behaviour ( p2| ) is already valid for g > 3, 
the MC results allow to determine the constant in (|2^), 
giving c = 80 from a two point fit. One should note that 
a prefactor of this magnitude has also been obtained for 
the low-temperature-behaviour of the correlation length 
of the nonlinear cr-modcl It is difficult to estimate, 
however, whether the asymptotic behaviour has already 
been reached at these values of g. Unfortunately the 
exponential increase of the correlation length does not 
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allow us at present to verify numerically our analytical 
result ( |22| ) for the g-dependence of the prefactor []l6| . In- 
deed the situation is again analogous to the much stud- 
ied 0(n) nonlinear a- model in d — 2. While the purely 
numerical constant equivalent to c(oo) in ( pl| ) was deter- 
mined approximately by Shenker and Tobochnik |i~8| ] via 
a Monte-Carlo RG method, the exact coupling constant 
dependence of the prefactor which follows from the two 
loop RNG has only been verified in recent years by 
extremely extensive numerical computations p9| . Since 
the range of physical interest in the single electron box 
problem is restricted to g values smaller than about five 
(i.e. R t > 1.3kf2, see Fig. 1) it is evident that our present 
numerical results fully cover the experimentally accessi- 
ble regime. 

In conclusion we have calculated the zero temperature 
smearing of the Coulomb staircase in the single electron 
box for arbitrary values of the tunnel conductance g. In 
contrast to previous work on this problem, which has con- 
centrated on the behaviour near n x = 1/2 and g <C 1 
we have discussed the slope at the center of the plateaus. 
It has been shown that this is determined by the he- 
licity modulus of an inverse square XY-model. Quanti- 
tative MC simulations in the physically relevant regime 
g = 0.1 .. .5 yield the expected crossover between per- 
turbation theory and the asymptotic behaviour. More- 
over, the two loop RNG calculation uniquely determines 
the analytical behaviour at large conductance and shows 
that previous instanton calculations are problematic. In 
the existing experiments || the conductance had a fixed 
value of order g exp rs 0.02, i.e. well in the perturbative 
regime. In fact the first oder correction in (10) was used 
to calibrate the slope at n x = to zero. To verify our 
results one therefore needs measurements with different 
and considerably larger values of g at temperatures where 
the thermal broadening is negligible. 
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FIG. 1. MC results for the zero temperature slope x of the 
Coulomb staircase at n x = as a function of the dimensionless 
conductance g. The dashed line is the perturbative result 
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Figure 1 
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